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THEORETICAL MODEL OF THE EFFECT OF CRACK TIP BLUNTING 
ON THE ULTIMATE TENSILE STRENGTH OF WELDS IN 2219-T87 ALUMINUM 


By 


Robert J. Beil 

Associate Professor of Engineering Mechanics 
Vanderbilt University 
Nashville, Tennessee 


ABSTRACT 


A theoretical model representing blunting of a crack tip radius 
through diffusion of vacancies is presented. The model serves as the 
basis for a computer program which calculates changes, due to successive 
weld heat passes, in the ultimate tensile strength of 2219-T87 aluminum. 
In order for the model to yield changes of the same order in the ultimate 
tensile strength as that observed experimentally, a crack tip radius of 
the order of .001 microns is required. Such sharp cracks could arise 
in the fusion zone of a weld from shrinkage cavities or decohered phase 
boundaries between dendrites and the eutectic phase, or, possibly, from 
plastic deformation due to thermal stresses encountered during the weld- 
ing process. 

Microstructural observations up to X2000 (resolution of about .1 

ed significantly under the influence of a heat pass, with the exception 
of possible small changes in the configuration of the interdendritic 
eutectic and in porosity build-up in the remelt zone. 
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INTRODUCTION 


The ultimate strength for a variable polarity plasma arc butt weld 
in 3/8- inch plate of 2219-T87 aluminum has been found! , by experiment, 
to increase when additional weld heat passes at lower power are made over 
the weld. Weld heat passes follow the same path as the original weld, 
but re-melt only a portion of the cast material in the fusion zone as 
shown in Figure la, b, c, d. Yield stress has been found to decrease 
during the same experimental procedure. The magnitude of these changes, 
obtained from tensile tests, is shown in Figure 2 !. 

This paper presents a theoretical model which indicates that, in the 
critical region where fracture tends to propagate, coarsening of crack tip 
radii by a diffusion process which takes place during weld heat passes is 
feasible and can contribute to increase in ultimate tensile strength. 


Characteristics of the physical process of fracture, especially in 
aluminum, have been studied, using the kinetic theory of materials, for 
over twenty years. Two review articles describe work at the A.F. Ioffe 
Physico-Technical Institute USSR Academy of Sciences, Leningrad, where over 
one hundred materials were tested, by 1970, for confirmation of the appli- 
cability of kinetic theory to the failure process. Materials tested in- 
cluded both metal and non-metal lie materials, single and polycrystalline 
materials, alloys, composites, and polymers. By 1980, there, the fracture 
process during creep tensile tests for aluminum had, through the use of 
relative density measurements accurate to 10 "^^ through the use of electron 
microscopy, x-ray diffraction, et. al . , been shown to include, at fracture, 
large numbers of sub-microscopic discontinuities (>10'/cm^) of length less 
than one-half micron and large numbers or microscopic discontinuities 
(>10°/cm!) of length greater than 1/2 micron in the lateral surface of the 
test specimen^. Creep tests^, were performed on polycrystalline aluminum 
(99.96%) in the temperature range 18-300°C and stress range 1-7 kg/mm^. 
Specimens were loaded for different times and then unloaded to study the 
build-up kinetics of microdiscontinuities. Researchers at AFIPTI report 
that rupture during creep develops from "atomic" cracks which arise after 
the application of the load, grow rapidly, become blunted, and may even 
change into stable void-shaped discontinuities not larger than !>«<«iby the 
end of the first stage. The mechanism for crack blunting, the mechanism 
for variations in microcrack concentration in the volume and in surface 
layers^, and the mechanism for coalescence of cracks and voids is a 


diffusion process most probably involving the emission of vacancies from 
breaks in continuity to sinks^. This paper restricts discussion to the 
role that diffusion of vacancies might have in blunting of crack tip radii. 
The reader is referred to related papers which discuss the role of 


dislocations, the role of original microporosity, and the role of hydro- 
static pressure in fracture kinetics which also contain materials which 


supports the feasibility of the process considered in this paper. 
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Figure la. Initial Weld Cross- \ 
Section. Current= 190 amps., 

Voltage= 31 volts. Weld Speed=7.3ipm 
1/16 dia. 2319 wire feed= 60 ipm/, 
Plate Thickness= 3/8 in. 


Figure lb. First Weld Heat 
Cross-Section. 

Current = 175 amps. 

Voltage = 28 volts 
Weld Pass Speed = 9 ipm. 

No Wire Feed 


Pass 
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Figure Ic. Second Weld Heat Pass 
Cross-Section . 

Current= 175 amps. 

Vo1tage= 28 volts 
Weld Pass Speed= 9 i pm . 

No Wire Feed 


Figure Id. Third Weld Heat Pass 
Cross-Section. 

Current= 175 amps. 

Voltage= 28 volts 
Weld Pass Speed= 9 i pm . 

No Wire Feed 





YIELD STRENGTH (ksi) ULTIMATE TENSILE STRENGTH (ksi) 


45 




Figure 2. Effect of Mutiple Heat Passes on the 

Strength of a Variable Polarity Plasma 
Arc Butt Weld in 3/8-inch 2219-T87 
Aluminum Plate. 
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PERCENT ELONGATION 
(2-inch gage length) 


THEORETICAL MODEL 


If the local temperature is high enough at a crack tip where the 
radius of curvature is small and the vacancy concentration is high, 
vacancies with diffuse to volumes of lower vacancy concentration. They 
may migrate through the bulk material to other voids or to flat surfaces. 
It is assumed that the elliptical crack considered is fairly flat at the 
ends of its minor axis so that vacancies will not need to travel further 
than half the crack length. The diffusion process increases the radius 
of the crack tip (blunting or coarsening). It is assumed that the 
vacancy flow will have little effect on the relatively flat portion of 
the crack. The diffusion concept is used in Appendix A to develop the 
differential equation 

- 


"*) fig dt 


where <1^ is the local ultimate stress 

o( is the coarsening angle firmed by tangents drawn 
to the assumed circular crack tip 

Co the equilibrium concentration of vacancies at 
the local temperature 

T is the local temperature 
is the crack dip radius 

Pj is the self diffusion coefficient of aluminum 

t is the time 

/U is the atomic volume of fee aluminum 

U is half the crack length 

V is the surface energy of fee aluminum. 


The above 
Appendix A, 


formula, the following differential 

dq^, „ dR 

~ z R ' 


equation, also derived in 


and a record of the history of the local temperature during the diffusion 
process^® form a system of equations which through iterative calculations 
determine the local ultimate strength and crack tip radius as a function 
of time. 
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CALCULATIONS 


Estimate of 

If a line is drawn tangent 
to the ellipse shown in 
Figure 3, then 



FIGURE 3 

Assuming that b ts y* and that L is much larger than the radius, R, of 
curvature of the ellipse at the point (L,0), so that 

L « X* and R = b, then 
TAN«t/2 = R/L 

Calculation of Atomic Volume 

The lattice parameter for fee aluminum is 4.04(10"®)cm. Hence, the 
cell volume is 65.94(10"^^)cm’ which accommodates four atoms. Thus the 
atomic volume is 16.48(10“24)cm3/atom. 

The Self- Diffusion Coefficient , Surface Energy , and Equilibrium Vacancy 
Concentration 

Volin and Balluffi^^ determined, as a consequence of study of the an- 
nealing kinetics of voids in aluminum, that the self diffusion coefficient 
of aluminum was 

Ds = .176 Exp (-1.31 ev/kT) cm^ 3"^ 

in the temperature range 85-209°C, which satisfactorily agreed with the 
work of Lundy and Murdock, and which agrees satisfactorily with the work 
by Bass describing the formation and motion energies of vacancies in 
aluminuml^. The surface energy of aluminum is taken as 1500 ergs/cm^ as 
in the paper by Volin and Balluffi even though calculations suggest that 
the figure might be slightly lower. 

Bass lists the formation energy of a vacancy in aluminum as .76 + 0.02 eV 
and the equilibrium vacancy in aluminum as 

C = Co EXP (.76 + .02 eV/kT) 
where Cq = 10^^ /16.48 cm“^ = 6.066 (10^^) cm“^ 

k = Boltzmann's constant = 8.611 (15“5) eV/°K 

and T is the temperature in degrees Kelvin. 
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Temperature History 


The Portion of the computer program which calculates the local tem- 
perature as a function of time was developed by Dr. A. C. Nunes, Jr, 
using a moving point source along with a moving line source to model the 
heat input. Non-conductive heat losses can be adjusted so that the width 
of the remelt zone of the fusion zone can be matched. 

Five points equally spaced along a line from the edge of the crown 
to the edge of the root on the opposite side of the fusion zone, as shown 
in Figure 4, were chosen as local points for calculation. Tensile test 
specimens show that fracture tends to occur along that line. Point A 
was taken to be just outside the remelt zone. The computer model for 
calculating temperature predicted temperatures too high at Point B when 
non-conductive losses were used to match the width of the remelt zones 
shown In Figure 1. Hence, higher non-conductive losses were assumed so 
that the calculated width of the remelt zone was reduced. For the newly 
calculated remelt zone, the circular part of the remelt zone was slightly 
larger, but still circular, and the shortest distance of the remelt zone 
from Point B remained the same. Temperatures for Points B, C, D, and E 
were calculated using the adjusted parameters. It is felt that this 
correction gave realistic results for the temperature distribution. 
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PERCENT INCREASE IN ULTIMATE 




WELD LOCATION o 
Crack Tip Radius= 20 A (,002^m) 

0 

Crack Length= 1000 A ( .100// m) 

Figure 4. Calculated Effect of Heat Pass on Local 
Ultimate Tensile Strengths Within the 
Primary Fusion Zone of a Variable Polarity 
Plasma Arc Butt Weld in 3/8-inch 2219-T87 
Aluminum Plate 
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DISCUSSION 


The amount of coarsening of the crack tip radius and the resulting 
change in ultimate strength during a weld heat pass is particularly sen- 
sitive to 

* the sharpness of the crack tip and, 

* the shape and magnitude of the local temperature 
versus time curve. 

Significant coarsening occurs, due to short heating times, only above 
temperatures of approximately 650°K or 680°K for crack tip radii of lo8 or 
50a respectively. Hence, significant improvement in local ultimate tensile 
strength takes place at only Points A, B, and possibly C of Figure 4. 

The local ultimate tensile strength of Points A, B, C, D and E are 
not likely to be identical immediately after the initial weld since the 
cooling rate is different for each of those points. This non-homogeneity 
does not complicate calculations since the configuration of the system of 
equations to be used and the iterative process used produce the relative 
change in ultimate tensile strength. Figure 4 illustrates the percentage 
relative change at the various points for a crack tip of 20a, about four 
times the lattice dimension of fee aluminum. The relationship of these 
changes to changes in ultimate strength observed experimentally is shown 
in Figure 5, where only the increase in the local ultimate strength at 
Point A matched the observed values, The 20 a initial crack tip radius, 
which provides significant increase in local ultimate tensile strength, 
is of a size large enough to be observed by the transmission microscope. 
Preliminary, but extensive, scanning electron microscopy was utilized at 
powers up to 2000X with good resolution to observe the microstructure of 
the weld cross section. Micrographs representing the general areas near 
Points A, B, C, D and E are shown in Figures 6a, b, c, d, and e for a 
cross section of the fusion zone on a specimen not subjected to weld heat 
passes. Obviously, 2000X is not enough magnification to see the thin 
sharp cracks pertinent to this discussion. Replication techniques will 
permit use of the transmission microscope, which can yield the desired 
magnification. The micrographs do reveal the dendritic structure in the 
cast material and the interdentritic eutectic which is a mixture of CuAl 2 
and <A. -aluminum. Microprobe analysis also reveals that a copper depleted 
zone occurs next to the eutectic similar to those mentioned by Doig and 
Edington.l^ One can realistically surmise that at the edge of the crown 
of the weld, where Point A is taken, and where the interdendritic eutectic 
appears normal to the surface of test specimen, see Figure 6a, cracks of 
submicroscopic dimensions may initiate in the copper depleted zone after 
welding. Coarsening of the crack tip radii of these types of cracks, large 
numbers of which may form during plastic flow prior to fracture under 
tensile loading, may contribute to the increase in ultimate tensile strength. 
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It appears, however, that the increase in ultimate strength may be 
due to more than one mechanism. Ultimate tensile strength increases 
whether or not the remelt zone overlaps the crown "corner” of the weld or 
previous weld pass. Hence, refinement of the dendritic structure in the 
remelt zone, due to more rapid cooling because of lower power, may also 
contribute to increase in ultimate strength. Further, reconfiguration of 
the interdendritic eutectic, possibly diffusing into the copper depleted 
zone or reforming through self diffusion, a slight variation noted in the 
micrograph study, may play a role in the process. 

Coarsening of crack tap radii, then, represents a feasible, but still 
not certain, or partial explanation, of multiple heat pass strengthening. 
The process Itself is of sufficient practical importance, in fracture 
mechanics, to be worth investiagation in its own right. Although progress 
has been made, further work is still necessary to pin down the mechanism 
of multiple heat pass strengthening. 
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PERCENT INCREASE IN ULTIMATE TENSILE STRENGTH 



Aluminum Plate. 
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Fioure 6a. Point A, Initial Weld. 
200X 

Alumina Polish 
Keller's Etch 


Figure 6b. Point B, Initial Weld. 
200X 

Alumina Polish 
Kel 1 er ' s Etch 
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Figure 6c. Point C, Initial Weld 
200X 

Alumina Polish 
Keller's Etch 


Fiqure 6d. Point D, Initial Weld. 
200X 

Alumina Polish 
Keller's Etch 
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Figure 6e. Po.int E, Initial Weld. 
200X 

Alumina Polish 
Keller's Etch 


RECOMMENDATIONS 


The following proposed theoretical and experimental research will 
help in identifying those processes which characterize fracture and in- 
crease ultimate strength in welded 2219-T87 aluminum. 

*..,to better understand the fracture process. 

Proposed experimental research should include 
microscopic examination, both low power, SEM, 
and replication TEM, of cross sections of the 
weld of tensile tested specimens. Test speci- 
mens should be loaded to a particular stress, 
held at that stress for a period of time, 
unloaded, and examined. This test procedure 
may allow identification of critical areas 
where cracks nucleate and characterize the 
growth pattern to failure. 

.establish a more accurate method to deter- 
mine the local temperature history. Higher 
order multipoles could improve the present 
technique. Finite element methods, allowing 
a variable thermal conductivity throughout 
the material, have been formulated for solution 
of boundary value and initial value problems 
to determine local temperature history such as 
needed in the paper. Since microstructural 
processes are so sensitive to local temperature 
history, computer software should be obtained, 
if it exists, or developed to determine local 
temperatures more accurately. Some measure- 
ments needed just to set up analyses such as 
loss parameters are not known. However, checks 
against the program could be verification by 
thermocouple readings near the weld and the 
pattern of the remelt zone after weld heat passes. 

*...to understand the effect of a refined dendritic 
pattern in the remelt zone after weld heat passes. 

A theoretical and experimental study should be 
made of the relationship of the size of the den- 
dritic structure with respect to ultimate stress. 

The size of the dendritic structure can be 
varied by varying the weld parameters of weld 
heat passes or by varying the quenching rate of 
the material during welding or weld heat passes. 
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I 


*...to better understand the role of the inter- 
dendritic eutectic. The diffusion coefficient 
of copper In aluminum Is known as a function of 
temperature. Hence, various, mechanisms, 
possibly self-diffusion within the eutectic, or 
diffusion of copper back Into the copper depleted 
zone as a function of time and temperature might 
be checked experlmentaly by heat treatment of cast 
material from the weld fusion zone and subsequent 
microscopic examination. 
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APPENDIX A 


EFFECT OF DEFECT EDGE RADIUS COARSENING ON ULTIMATE TENSILE STRENGTH 

The preliminary version of the analysis was prepared by 
Dr. A. C. Nunes, Jr, 


Ductile fractures generally occur by coalescence of voids within the 
fracturing material . Voids often form around defects such as cracks or 
weakly bonded second phase particles. The limitation on the ultimate ten- 
sile stress, is set by the critical stress for propagation of internal 
voids. The maximum stress supportable by the defects within the metal 
should be that for initiation of a void from the defect. 


If the defect has a sharp edge of suitable orientation, the concentra- 
tion of stress at the edge should be where the void starts. The stress 
concentration for an ellipsoidal crack of length iLand radius of curvature 
^ at the crack edge raises the nominal stress <r^to a higher local effec- 
tive stress according to the Inglis relation: 



Assuming that the nominal stress becomes the ultimate tensile strength 


at a critical value of cl 




then 


r ^ V? 

U X ^ 


If the radius of the sharp edge of the defect changes, the ultimate 
tensile strength changes, then, according to the relation: 

d <T^ _ Cu, 
df> ’ 

If two faces of a defect come together at angle <x. to a radius a as 
shown in Figure A-1, then a loss of defect volume at the edge would result 
in a change of radius Ap . 



FIGURE A-1. EDGE OF A DEFECT SHOWING CHANGE OF RADIUS 
OCCURRING WITH LOSS OF VOLUME AV 
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The shaded area, can be shown to be 

{2P^e 

so that ^ * < 2 -^^ {<«,r (‘‘/jl) -C iT-ot^/fc} 

disregarding second order terms. 

The kind of volume reduction at the defect edge shown in Figure A-1 
takes place because of a difference in chemical potential, , for solutes 
(or vacancies) on surfaces of different curvature. When a solute atom or 
vacancy moves from a cylindrical surface of radius to one of radius 
causing a volume loss Jl. to the first cylinder and a^volume gain of Ti- 
to the second, a net change A^in surface area results 

A A - -n. C 

If the surface energy of the cylinders is if' per unit area then the 
chemical potential difference between the surfaces is: 

A/4 * yTL(-^--}r') 

Assuming that the activity of the solute or vacancy is proportional 
to its concentration, then the equilibrium concentration Cf at the curved 
surface is related to the equilibrium concentration c^at a flat surface 
according to the equation; 

It is further assumed here that the surface reaction rates are fast 
enough so that equilibrium concentrations are maintained at the surfaces 
and that the relatively lengthy time required for diffusion of solute atom 
or vacancy from regions of high concentration to regions of low concentration 
is what holds back the coarsening process. 

Pick's first law relating mass flux ^ in atoms or vacancies per area 
per unit time to the concentration gradient in the radial direction r is 

assumed to hold: : .n^ 

a - ar 

where 0 is the solute or vacancy diffusivity. 

A steady cyl indrical ly symmetrical flow of /w atoms or vacancies per 
unit time yields a flux; /it 

^ ' jCtrrTL 


Which, with Pick's first law, requires a concentration Cr varying with the 
radius according to the relation: ^ ^ 

Cf.~ " aitA D f 


Where is the concentration at radius /» . 
Hence : 

1 
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Assuming that there is so much more flat surface ( than curved 
edge that the bulk mean solute or vacancy concentrate Cois determined by 
the flat surface, then: 


for large r and 




Cp-Cr ^ 

c 

The loss of volume v from the sharp edge of a defect is given by; 

^ = -A- (. IT- ^ ^ ^ji. 


o r 


Iv s jL oc^ ( 

d t '-ff') 

Where Lis half the crack length and assumed an approximate maximum upper 
bound distance to a flat surface. 

The ultimate tensile strength then varies according to the relation: 

^ 

^ dfc'clf ’dv'dt 

so that 

-u D 


d 




6t 
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I 


ro 


cn 


I REM tit RhPID DIFFUSSIOH ANNEALING tt* 

I PEN itt INCREASE IN ULTIMATE STRENt:TH ttt 
I CLEhR 

I ^ PRINTED OUTPUT WANTED > 0 OTHERWISE." 

PRINT "THIS PROGRAM CALCULATES INCREASED ULTIMATE" 

PRINT "STRENGTH DUE TO RAPID DIFFUSION ANNEALING" 

PRINT BLUNTINij) OF MICROCRACKS. TEMPERATURES AT" 

FRUIT "A MATERIAL POINT APE TO BE ENTERED FOP DISCRETE" 

PRINT CALCULATED FROM THE ATTACHED SUBROUTINE." 

PRUa ’•INPUT 0 TO SKIP ALL PRINTING EXCEPT THE LAST " 

PRINT "INPUT I TO GET TEMP > POINT DATA. ET AL " 
lUPiJJ ZA 

PFINT "ENTER 0 IF THE SUBROUTINE IS TO BE CALLED " 

PRINT "ENTER 1 OTHERWISE" 

INPUT B 

PRINT "If^PUT THE NUMBER OF HEAT PASSES NOT INCLUDING THE WELD PF 
INPUT WP 

PRINT "INPUT NPTS-l ACROSS THE DIAGONAL " 

IfiPUT LU 

DELETE UP.RP.TU.TR 
DIM UP^'iJP.LU >.RP<WP.UO 
D I M TU* WP , LU ;• .. TR< WP . LU > 

FOR 11=0 TO WP 
IF B=1 THEfI 2030 


200 REM i^t Hh2 TEMPERATURE HISTORY ttt 
210 REM tit data If^PUT ttt 

230 PRINT "PLEASE SELECT THE NUMBER OF THE METAL TO 

240 PRINT "FROM THE LIST BELOW . "vPRINTvPRINTKPRINT 

250 PRINT "0 UNLI-STEO METAL"\PRINT 

2C0 PP.UiT "1 2219 ALUM I HUM ‘'SPRINT 

2Z0 PRIfU "2 302 stainless STEEL"\PRINT 

230 PRINT "3 321 STAINLESS STEEL"\PRINT 


BE WELDED 


APPENDIX B 


290 PRINT "4. INCONEL 718"' PRINT 

300 INPUT Q 

310 IF Q=1 THEN 550 

320 IF 0=2 THEN 570 

3-^0 IF 0=3 THEN 590 

340 IF 0=4 THEN CIO 

350 PRINT "WHhT IS THE METAL TO BE WELDED?" PRINT' PR INT- PRINT 
3C0 INPUT Mi- 

370 PRINT "PLEASE WRITE IN THE LOWER AND UPPER LIMITS OF THE MELTING" 
3S0 PRINT "TEMPERATURE RANGE OF THE ALLOY TO BE WELDED. USE UNITS OF" 
3*^0 PRINT "DEGREES FAHRENHEIT SEPARATE THE TWO MALUES BY A COMMA." 

400 PRINT "NOTE^ A GOOD SOURCE REFERENCE FOR THIS KIND OF DATA IS THE" 
410 PRINT "AEROSPACE STRUCTURAL MATERIALS HANDBOOK . " ‘.RRINT'^PRINT -PRIMT 
420 INPUT LT.UT 

430 PRINT "WHAT IS THE THERMAL CONDUCTIUITY OF THE ALLOY IN BTU’S PER" 
440 PRINT "HOUR PEP FOOT PER DEGREE FAHRENHEIT?" 'PR I NT 'PR I HT\PR I NT 
450 INPUT Kl 
4C0 LET K1=K1T2 441E-05 

470 PRINT "WHAT IS THE SPECIFIC HEAT OF THE ALLOY IN BTU’S PER LB" 

430 PRINT "PER DEGREE FAHRENHEIT?"' PRINTvpRINT' PRINT 
490 INPUT Cl 

500 PRINT "WHAT IS THE DENSITY OF THE ALLOY IN LBS PER CUBIC INCH?" 

510 PR I NTvpR I NT' PR I NT 
520 INPUT PO 

530 LET A1=K,1^R0^C1T56.83 
540 GOTO 630 

550 LET UT=1190'LET LT=1010''LET Kl=2 . 44E-03'LET Al=5 56 
560 LET Mi ="22 19 ALUMINlJM"\GOTO 630 

570 LET UT=2650 LET LT=2550'LET fa =4 . 27E-04 'LET Al= 5 

530 LET Mf="302 STAUfLESS STEEL" GOTO 630 

590 LET UT=2550' LET LT=25O0 LET Kl=4 O3E-04 -LET Hl= 66 

600 LET Mi="321 STAINLESS STEEL"'GOTO 630 

610 LET UT=243? LET LT=230O LET Kl=3 .83E-04'.LET Al= 71 

620 LET M$="IMCOtfEL 718" GOTO 630 

630 PRINT "THE METAL TO BE WELDED IS NOW CHARACTERIZED WE NOW" 



640 PRINT "TURN TO THE WELDING PROCESS PARAMETERS. “nPR I MT\PR I NT vPRINT 
650 PRINT "WHAT IS THE AMBIENT TEMPERATURE OF THE METAL IN DEGREES” 

660 PRINT "FAHRENHEIT? IF THE METAL IS PREHEATED. GIUE THE PREHEAT” 

670 PRINT "TEMPERATURE AS THE AMBIENT TEMPERATURE . ”nPR1 NT nPR I MT^PR I NT 

6S0 INPUT TO 

760 PRINT "PLEASE SPECIFY WELD POWER IN KILOWATTS.” 

770 PRIMT\PRINT'.pRINT 
780 INPUT PO 

780 PRINT "PLEASE SPECIFY WELD SPEED IN INCHES PER MINUTE.” 

800 PR I NT' PR 1 NT 'PR I NT 
810 INPUT U1 

^ S20 PRINT "WHAT PERCENT OF TOTAL BEAM POWER IS LOST FROM THE WELD” 

S -30 PRINT "PUDDLE DUE TO PROCESSES OTHER THA^^ CONDUCTION BY THE PLATE?" 

►-I 840 PRINT “THESE LOSSES INCLUDE RADIATION AND METAL EUAPORATION FROM” 

850 PRINT "THE UICINITY OF THE PUDDLE . " sRRIHTnPRIHT'PRINT 

860 INPUT FI 

861 PRINT "INPUT THE NEW NON-CONDUCT I UE PERCENT LOSSES FOR DEEP PTS.“ 

862 INPUT NC 

863 CC=F1 

370 PRINT "WHAT PERCENT OF THE REMAINING BEAM POWER IS ABSORBED CLOSE” 

880 PRINT "TO THE METAL SURFACE SO AS TO FORM THE EB WELD NAILHEAO?” 

890 PRINT "THIS WOULD BE THE PERCENTAGE OF THE BEAM CURRENT LACKING” 

900 PRINT "SUFFICIENT POWER DENSITY TO UAP0RI2E THE METAL.” 

910 PRINT' PR I NT' PR I NT 
920 INPUT F2 

940 DIM KDri50).TSCl50)JJS-r. 150).R<150) 

950 IF II^;>0 THEN 1140 

960 PRINT "PLEASE SELECT LENGTH OBSEfUED 0.5. 2.5. 5. 10. 25 INCHES.” 

970 PRINT 'PRINT' PRINT 

930 INPUT XB LET XB=XB/5nLET XA=-4*XB 

990 PRINT “PLEASE SELECT MAXIMfJM TEMPERATURE^ 1000. 2000. 5000. 10000 DEG 
F.” 

1000 PRIHT\RRINT\PRIHT 
1010 INPUT TG 

1020 PRINT “INPUT HALF THE CRACK LENGTH IN CENTIMETERS.” 


1030 INPUT L 

1040 PRINT "INPUT THE LOWER MELTING TEMRERATUREC DEG . KELUIN) " 

1050 INPUT LQ 

1060 PRINT "INPUT THE RADIUS OF THE TIP OF THE CRACK<CM>.” 

1070 INPUT R<0> 

1030 REM R IS ASSJJMED MUCH LESS THAN L 
1090 print "INPUT THE ULTIMATE STRESS.” 

1100 INPUT uS':o:> 

1102 FOR MM=0 TO LU 

1103 UPC 0. MM >=U3<0> 

1104 PP<0.MM>=R<0> 

1105 NEXT MM 

1106 PRINT'PRINT 

1110 PRira "INPUT ONE HALF THE ROOT WIDTH. OJC HALF THE CROWN WIDTH." 
1112 PRINT "AND THE PLATE TH I CKNE3SC DIMENSIONS IN INCHES)." 

1120 INPUT BH.TH.Wl 

1125 DELETE VY. 2Z. MU.MR. SU, SR 

1130 DIM YY< LU > . ZZ<. LU > . MU< LU ) . MR< LU > . SUC LU ) . SRC LU ) 

1140 FOR MM=0 TO LU 

1145 IF IIO0 THEU U3C O >=UP(I I . MM > 

1146 IF IK >0 THEN R< 0 )=RP( 1 1 . MM > 

1150 Y=-BH+< TH^BH >TMM/LU 

1160 Z=-WUC Y-THVCBH+TH> 

1170 YYCMM>=Y 
11 SO Z2<' MM)=Z 
1190 GOSUB 3000 
1200 whit 1500 
1210 PAGE 
2020 GOTO 2130 

2030 PRINT "INPUT ONE LESS THAN THE NUMBER OF DISCRETE" 

2040 PRIf^T "TEMPERATURES TO BE CONSIDERED." 

2050 INPUT M 

2060 DIM KDCrD.TSCM) 

2070 DIM IJS<M>.R<M> 

2030 FOR 1=0 TO f1 



PRINT "irFUT TE^FER^^TlJRE<OEG. KELUIN >>TIME< SECONDS)” 

^iOO INPUT KOa >..TS-; I > 

21 lu fiE:.T I 

2120 GOTO 2220 

2170 FOP I-O TO 150 

2140 Th-: I ..■=5rTh4< I V9+255.77 

2150 me: :t 1 

21 GO DT=< 2T3>AU 
2170 M=150 
2100 FOR 1=0 TO 150 
< 2190 KO*: I .'=Th( 150-1 > 

2200 ne; :t i 

^ 22S0 hL=2 # ATN< R'.; O >y'L > 

, 2-300 PI=4t.HTfKl > 

ro 2310 0M= 1 6 . 434” 1 0 • ( -24 ) 

2740 page 

2741 IF PU=0 THEN 2420 

2350 PRINT ”Z=”;2; ”V=”;Y. ”L=”;L>‘'R<0>-”;P<0> 

2760 PRINT ”AL=*’.AL;”U3<0>=";US<0) 

2770 PRINT 

2780 PRINT "KDCD DS RA RR US< I > RCI )•' 

2390 PRINT "A CO S D LG” 

2400 PRINT 

2410 PRINT 

2420 FOR 1 = 1 TO ri 

2470 IF KOC I >^=LQ THEN 2470 

2440 US'- I '^=UR'' 0 , Mfl > 

2450 P*: I :-=PP^O.riM> 

2452 IF 2A=0 THEN 2460 

2455 PRINT 1 1 , I / Mf1;K0< I R< I-l USU-1 > 

2460 GOTO 2610 

2470 IF B=0 THEN LET DS=DT 

2480 IF B=1 THEN LET DSaTS< I >-TS< I-l ) 

2490 A=-' Df1*.< PI-AL >TOS >/< 44:< SIN< < PI-AL V2 )AOSC < PI-AL V2 ) >-< < PI-AL ^/2 ^ > 
2500 C0=6 O66410-< 22 )4^EKP< -3473/KD< I > ) 


2502 IF ZA=0 THEN 2510 

2505 PRINT II : I MH,KD- I :>..Pa-l :v.U3< I-l > 

2510 S=E: :P< 15Cu:ifOH. Y l 7S1*:10 ‘'-16 >>.' P-^ I-l VKOr I >> 
2-520 D= 1 76 ♦ E; IP'. -15213- rO< I > ) 

2570 LG=LOG'-' L- 'P«.' I - 1 > .) 

2540 RH=H»C0tS4D.'':LG4R<; I-l >'2> 

2550 RR=24PA 

2560 R< I >=R*' I-l >+RP4P<: 1-0 
2570 US^. I >=US< I-l >+PA4U3< I-l ) 

2575 IF PiJ=0 THEN 2670 
2580 PRINT 

2590 PRINT KD< I >. DS; PA; RR; US< I y; R< I ) 

2600 PRINT h;C0;S>D;LG 
2610 PRINT 

2615 IF PU=0 THEN 2630 
2620 PRINT FDCD.LO.I 
2670 NEXT I 

2632 IF II=NP THEN 2640 

2635 UPC 1 1 + 1 . firi >=IJS< 1 50 

2636 RPC 1 1 + 1 . riri >=RC 1 50 > 

2640 Mucm:>=cu3c i50>-uP''0.riri>A4jp<0,riM> 

2650 riPC MM >=C PC 1 05 >-PpC O . MM VRPC 0 . MM > 

2652 .SUCf1f1)=nUCMM> 

2654 SRC MM )=MRC TIM > 

2655 TLIC 1 1 . MM >=MUC MM > 

2656 TRC 1 1 > MM >=MPC MM > 

2660 NEXT MM 

2665 WAIT 1500 

2700 SU=C SIJ-^ MIIK SU > > >/C MA'XC SU )-MINC SU > > 

2710 SR=C SR-C MIHC SR > > VC MAXC SR VMINC SR > > 

2720 PAGE 

2724 IF II=WP THEN 2730 

2725 IF 2m=0 then 2300 

2730 PRINT "MAX DUT3/UTS= ”>MA>:CMU> 

2740 PRINT ”MIN DUTS/UTS= ”.MINCf1U> 




] 

ro 

00 


2750 WINDOW -BH/TH.0.1 
2760 UIEWPORT 250.750/ 100^ 6ei0 
2770 3ETGR TICS 5/5/GPrtT 5/5 
2780 XYPLQT YY/SU 
2750 PRINT "^C-W'-'C^L” 

2800 NEXT II 
280'^ PAGE 

2310 PRINT "DTS/UTS FOR WP=0/3” 

2812 PRINT ‘•E%*‘D’‘/"C"/“B“/’'A'* 

ilie PPIHT°TUf. Jm.Z, TU< 00, 1 >, TU< 00,2), TU( 00, 3 ), UK 00, 4 ) 
2818 PRINT 

2820 NEXT OQ ^ 

2340 PRINT "DR-'R FOR HP=®'3 

2842 PRINT "E" , "0”, "C", "B‘ , A 

Ittt PRIN?'’tR< 00 , Z, TR< 00, 1 ), TR< 00,2 >, TR< 00, 3 >, TR( 00, 4 > 

2843 PRINT 
2350 NEXT QQ 
2855 WhIT 1500 
2860 PRU4T 
2365 WAIT 1500 
2900 STOP 

§000 REM computation OF HA2 TEMPERATURE HISTORY *»* 

3005 IF II=GP THEN 30O9 

3006 IF f1M=LN THEN 3009 

3007 Fl-NC 

3008 GOTO 3010 

3009 F1=CC 

3010 Pl=P04^< 1-^1/100 
3015 PRINT II/riM/UU/Fl 
3020 P2=PU< 1“F2/100> 

3030 Pl=Pl-P2 

3040 LET N=150 


3050 DIM XX<N>/TA<M>/TEvN> 

3060 LET TE=LT 

3070 FOR 1=0 TO H 

3030 LET XXC I >=:<A+<XB-^:A >/HtI 

3090 NEXT I 

3110 FOR 1=0 TO N 

3120 LET X=XXd > 

3130 GuSUB 3290 
3140 LET TAd)=M2 

3145 IF TA<I)>LT THEN LET TA< I )=LT+10 
3150 NEXT I 

3160 PEH PRINTOUT OF RESULTS m 
3170 PAGE 

3175 IF 2A=0 THEN 3230 

3130 PRINT "HELD TEMPERATURE PROFILE "/ Y. INCHES FROM CENTERLINE 
3190 PRINT "PROFILES TAKEN AT " / Z/WUIOO; *' *; PLATE DEPTH." 

3200 FOR 1=0 TO 10 PRINT NEXT I 
3210 PRINT " TEMP"'PRIHT'PRINT 
3220 PRINT "<DEG F /' 

3230 FOR 1=0 TO 16' PRINT' NEXT I 

3240 PRINT " DISTANCE FROM HEAT SOURCE < INCHES >"/ 

3250 WINDOW XA/XB/O/TG 

3260 SETGR WIMDOW/TICS 5/ lO/ 5/ 5/ GRAT 6>6/3/3 
3270 XYPLOT XX/ TA 
3275 PRINT "-C W'C-L" 

3230 RETURN 

3290 PEM tt* SUBROUTINE TO COMPUTE TEMPERATURE OF HEAT SOURCE aPRhY 

3300 LET Z2=Z 

3310 LET M2=0 

3320 LET D3=0 

3330 LET D2=2»W1 

3340 GOSUe 3650 

3350 G03UB 3510 

3360 LET M2=M1+T0+T2 

3370 LET M3=M2 



I 

ro 


3330 LET D3=D3+D2 
7330 LO’ZUB 3430 

3400 IF <M3-M3><01*M2 THEN GOTO 3420 
3410 GOTO 3370 
3430 RETURri 

3430 PEM 44 ♦ SUBROUTINE TO ADD NEXT TWO HEAT SOURCES ttt 

3440 LET Z2=D3+Z 

3450 GOSUe 3b50 

3460 LET M2=M2+m 

3470 LET 22=D3“Z 

3480 GOSUB 3650 

34?0 LET M2=N2+M1 

3500 RETURN 

3520 PEft DISTRIBUTION 444^ COMPUTE LHC HEAT SOURCE TEMPERATURE 

3530 LET 10=0 

3540 LET T2=0 

3550 LET AR=Ul/2/A14X 

3560 LET Z3=Ul/2 'A14SQR<X4X+Y4y > 

3570 IF AR>88 THEN -3640 

3580 GOSUB 3730 

3580 IF AR<0 THEN 3629 

3600 LET T2=P2/Wl/6.2332/Kl/EXP'^AR)4K0 

3610 GOTO 3640 

3620 LET hP=“AR 

3630 LET T2=P2."W1/6.2332^'K14EXP< AR>4K0 
3640 RETURN 

3650 REM 444 SUBROUTINE TO COMPUTE TEMPERATURE 444 
3660 LET S1=SQR*'X4X-*'Y4V+Z2422:> 

3670 IF SI=0 THEf< LET Ml=UT4lOO0^IF 31=0 THEN 3720 
3680 LET hR=U1/2/A14<SH-X> 

3690 IF AP>88 THEN LET M1=0 

3700 IF AR>83 THEM GOTO 3720 

3710 LET Ml=Pl/2/3. 1416/KI/S1/£XP-^AR> 

3720 RETURN 


3730 REM 444 MODIFIED BESSEL FUNCTION. SECOND KINO. ZEROTH ORDER 444 

3740 REM 441 POLYNOMIhL hPPRO: 'IMATION 444 

3750 IF Z3;2 THEN 3820 

3760 GOSUB 3^70 

3770 LET G2=Z3423/4 

3780 LET K0=< < 7 . 4E-064G2+ 1 . 075E-G4 >4G2+2 62693E-03 >IG2 

3790 LET K0=- <K0+ . 0348859 >4G2+ 2306976 >4G2 

3300 LET K0=<<KO+ . 4227842 >4G2- 5772157 >~.54L0G(:G2 >410 

3310 GOTO 3860 

3320 LET G2=2/Z3 

3S30 LET K0=< < 5 . 3208E-'044G2-2 . 5154E-03 >4G2-^5 . 87872E-83 >4G2 

3340 LET K0=( < K0~ .01062446 >4G2+ . 02185563 >4G2 

3350 LET K0=< (K0- . 07832358 >4G2-H 253314 >/SOR< 23 >/EXP< Z3 > 

3860 RETURN 

3870 REM 444 MODIFIED BESSEL FUf4CTI0H. FIRST KIND. ZEROTH ORDER 444 
3330 REM 414 POLYMOMIAL APPROXIMATION 444 
3890 IF 23 >3. 75 THEN 3958 
3900 LET G 1=23423 

3910 LET I0=<<5 923979E-104G1+6 56017E-08>4G1*^6.80123E-O6..4G1 
3920 LET I0=C < 10+4 . 3394E-84 )4G1+ 0156252 >4G1 
3930 LET I0=< I0+.25>4G1+1 
3940 GOTO 3990 

3950 LET I0=(< 153.445/23-171 .822 >/23+73. 2919 V23 

3960 LET I0 =< (. 10-15 . 2595 >/23+l 81 193 >/23 

3970 LET I0=< <10- . 0830909 >/23+ 0316355 >/23 

3980 LET I0=< < 10+ . 0493222 >/23+ . 3989423 >4EXP< 23 >/SOR< 23 > 

3990 RETURN 

READY 

4 


I 


10 REM ttt RAPID DIFFUSSION ANNEALING ttt 
20 REM ttt INCREASE IN ULTIMATE STRENGTH tt* 

30 CLEAR 

35 PRINT "INPUT 1 IF MAX PRINTED OUTPUT WANTED. 0 OTHERWISE." 

36 INPUT PU 

40 PRINT "THIS PROGRAM CALCULATES INCREASED ULTIMATE" 

50 PRINT "STRENGTH DUE TO RAPID DIFFUSION ANNEALING" 

52 PRINT "< BLUNTING) OF MICROCRACKS. TEMPERATURES AT" 

54 PRINT "A MATERIAL POINT ARE TO BE ENTERED FOR DISCRETE" 

56 PRINT "TIMES OR CALCULATED FROM THE ATTACHED SUBROUTINE." 

60 PRINT 

62 PRINT "INPUT 0 TO SKIP ALL PRINTING EXCEPT THE LAST." 

63 PRINT "INPUT I TO GET TEMP . POINT DATA. ET. AL." 

64 INPUT 2A 

70 PRINT "ENTER 0 IF THE SUBROUTINE IS TO BE CALLED." 

72 PRINT "ENTER 1 OTHERWISE" 

80 INPUT B 

81 PRINT "INPUT THE NUMBER OF HEAT PASSES NOT INCLUDING THE WELD PASS." 

82 INPUT WP 

83 PRINT "INPUT NPTS-l ACROSS THE DIAGONAL." 

84 INPUT LU 

85 DELETE UP.RP.TU.TR 

86 DIM UP<P4P.LU>.RP<WP.LU) 

87 DIM TU<WP.LU).TR<WP.LU) 

88 FOR 11=0 TO WP 
90 IF B=l THEN 2030 

200 REM HAZ TEMPERATURE HISTORY ttt 
210 REM ttt DATA INPUT ttt 

230 PRINT "PLEASE SELECT THE NUMBER OF THE METAL TO BE WELDED" 

240 PRINT "FROM THE LIST BELOW . "^PRINT^PRINT^^PRINT 

250 PRINT "0. UNLISTED METAL"nPRINT 

260 PRINT "1. 2219 ALUM I HUM "\PR I NT 

270 PRINT "2. 302 STAINLESS STEEL"\PRINT 

280 PRINT "3. 321 STAINLESS STEEL"\PRINT 
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290 PRINT “4. INCONEL 718"nPRINT 
300 INPUT 0 . 

310 IF Q=1 THEN 550 

320 IF Q=2 THEN 570 

330 IF 0=3 THEN 590 

340 IF 0=4 THEN 610 

350 PRINT "NHAT IS THE METAL TO BE WELDED?“'PRINTnPRINTnPRINT 
360 INPUT Mt 

370 PRINT “PLEASE WRITE IN THE LOWER AND UPPER LIMITS OF THE MELTING" 
380 PRINT "TEMPERATURE RANGE OF THE ALLOY TO BE WELDED. USE UNITS OF" 
390 PRINT "DEGREES FAHRENHEIT. SEPARATE THE TWO 'JALUES BY A COMMA." 

400 PRINT "NOTE: a GOOD SOURCE REFERENCE FOR THIS KIND OF DATA IS THE" 
410 PRINT "AEROSPACE STRUCTURAL MATERIALS HANDBOOK. "xPRINT^PRINT^PRINT 
420 INPUT LT.UT 

430 PRINT "WHAT IS THE THERMAL CONDUCTIUITY OF THE ALLOY IN BTU'S PER" 
440 PRINT "HOUR PER FOOT PER DEGREE FAHRENHEIT?"nPRINT\PRINTnPRINT 
450 INPUT K1 
460 LET Kl=Kl:»:2.44lE-05 

470 PRINT "WHAT IS THE SPECIFIC HEAT OF THE ALLOY IN BTU'S PER LB" 

480 PRINT "PER DEGREE FAHRENHEIT?"\PRINTvPRINTnPRINT 
490 INPUT Cl 

500 PRINT "WHAT IS THE DENSITY OF THE ALLOY IN LBS PER CUBIC INCH?" 

510 PRINTxPRINTnPRINT 
520 INPUT RO 

530 LET A1=K1-^R0/Cl:lt56.83 
540 GOTO 630 

550 LET UT=1190vLET LT=1010 nLET K1=2.44E-03'4_ET Al=5.56 
560 LET M$="2219 ALUM I HUM "'-GOTO 630 

570 LET UT=2650 'LET LT=2550nLET K1*4.27E-04\LET A1= .5 

530 LET Mf="302 STAIhH-ESS STEEL"\GOTO 630 

590 LET UT=2550NLET LT=2500xLET K1*4.03E-04\LET Al=.66 

600 LET M*="321 STAINLESS STEEL"\GOTO 630 

610 LET UT=2437 -LET LT=2300vLET Kl*3.83E-04sLET Al=.71 

620 LET M$=" INCONEL 718"' GOTO 630 

630 PRINT "THE METAL TO BE WELDED IS NOW CHARACTERIZED. WE NOW" 
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640 PRINT "TURN TO THE WELDING PROCESS PARAMETERS. "nPRINTxPRINTxPRINT 
650 PRINT "WHAT IS THE AMBIENT TEMPERATURE OF THE METAL IN DEGREES" 

660 PRINT "FAHRENHEIT? IF THE METAL IS PREHEATED, GI»v»E THE PREHEAT" 

670 PRINT "TEMPERATURE AS THE AMBIENT TEMPERATURE. "xPR I NT%PR I NTvPR I NT 
6S0 INPUT T0 

760 PRINT "PLEASE SPECIFY WELD POWER IN KILOWATTS." 

770 PRINTnPRINTnPRINT 
780 INPUT P0 

790 PRINT "PLEASE SPECIFY WELD SPEED IN INCHES PER MINUTE." 

800 PRINTxpRIHTvPRIHT 
810 INPUT U1 

820 PRINT "WHAT PERCENT OF TOTAL BEAM POWER IS LOST FROM THE WELD" 

830 PRINT "PUDDLE DUE TO PROCESSES OTHER THAN CONDUCTION BY THE PLATE?" 
840 PRINT "THESE LOSSES INCLUDE RADIATION AND METAL EUAPORATION FROM" 

850 PRINT "THE UICINITY OF THE PUDDLE . "xPRINTsPRINTvPRINT 

860 INPUT FI 

861 PRINT "INPUT THE HEW NON-CONDUCT I UE PERCENT LOSSES FOR DEEP PTS." 

862 INPUT NC 

863 CC=F1 

870 PRINT "WHAT PERCENT OF THE REMAINING BEAM POWER IS ABSORBED CLOSE" 

880 PRINT "TO THE METAL SURFACE SO AS TO FORM THE EB WELD NAILHEAO?" 

890 PRINT "THIS WOULD BE THE PERCENTAGE OF THE BEAM CURRENT LACKING" 

900 PRINT "SUFFICIENT POWER DENSITY TO UAPORIZE THE METAL." 

910 PRINTvPRINT\PRINT 
920 INPUT F2 

940 DIM KD( 150 >, TS< 150 ), US< 150 R< 150 ) 

950 IF IIO0 THEN 1140 

960 PRINT "PLEASE SELECT LENGTH OBSER<viED= 0.5, 2.5, 5, 10, 25 INCHES." 

970 PRINT'.PRINT\PRINT 

980 INPUT XB' LET XB=XB/5\LET XA=-44cXB 

990 PRINT "PLEASE SELECT MAXIMUM TEMPERATURE = 1000, 2000, 5000, 10000 DEG 
F." 

1000 PRIHT^PRINT' PRINT 
1010 INPUT TG 

1020 PRINT "INPUT HALF THE CRACK LENGTH IN CENTIMETERS." 
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1030 INPUT L 

1040 PRINT "INPUT THE LOWER fCLTING TEMPERATURE< DEG . KELMIN>." 

1050 INPUT LQ 

1060 PRINT "INPUT THE RADIUS OF THE TIP OF THE CRACK<CM>." 

1070 INPUT R<0> 

1080 REM R IS ASSUMED MUCH LESS THAN L 
1090 PRINT “INPUT THE ULTIMATE STRESS.'* 

1100 INPUT US<0> 

1102 FOR MM=0 TO LU 

1103 UPC0.MM)*US<0.> 

1104 RP<0.MM>*R<0) 

1105 NEXT MM 

1106 PRINTnPRINT 

1110 PRINT “INPUT ONE HALF THE ROOT WIDTH. ONE HALF THE CROWN WIDTH. 
1112 PRINT “AND THE PLATE TH I CKNESS< DIMENSIONS IN INCHES)." 

1120 INPUT BH.TH.Wl 

1125 DELETE YV,22.MU..MR.SU.SR 

1130 DIM 'mLU)/ZZ<LU>.MU<LU).MR<LU).SU<LM).SR<LM) 

1140 FOR MM=0 TO LU 

1145 IF IIO0 THB< US<0)=UP< II.MM) 

1146 IF IIO0 THEN R<0>=RP< II.MM) 

1150 Y=-BH+< TH+BH >i|cMM/LU 

1160 2=-Wl3|c<Y*-TH>/<BH+TH> 

1170 YY<MM>=Y 
1180 Z2<MM)=Z 
1190 GOSrjB .3000 
1200 WAIT 1500 
1210 PAGE 
2020 GOTO 2130 

2030 PRINT "INPUT ONE LESS THAN THE NUMBER OF DISCRETE" 

2040 PRINT "TEMPERATURES TO BE CONSIDERED." 

2050 INPUT M 

2060 DIM KO<M).TS<M) 

2070 DIM US<M>/R<M> 

2080 FOR 1=0 TO M 



2090 PRINT "INPUT TEMPERATURE': DEG . KELUIH T I f1E< SECONDS)" 

21V30 INPUT KDa>..TS':;i) 

2110 NEKT I 

2120 GOTO 2220 

2130 FOR 1=0 TO 150 

2140 TR-: I )=5TTA< I >/9+255 . 37 

2150 ne:<;t I 

2160 DT=<2:tXB>/Ul 

2170 M=150 

2180 FOR 1=0 TO 150 

2190 KDC I ;>=Th< 150-1 > 

2200 NEXT I 

2290 AL=2*ATN<:.R(;0VL) 

2.300 PI=4#ATH<1) 

23 1 0 OM= 1 6 . 48* 1 0' 6 -24 

2340 PAGE 

2341 IF PU=0 THEN 2428 

2350 PRINT "2=";2;"Y=";Y;"L=";L^"R<0)=";R<0> 

2-360 PR I NT " AL= " ; AL ; " US< 0 >= " ; US< 0 > 

2370 PRINT 

2380 PRINT "KD(. I > DS RA RR US( I ) R< I )" 

2390 PRINT "A C0 S D LG" 

2400 PRINT 

2410 PRINT 

2420 FOR 1=1 TO M 

2430 IF KD<IX=LQ THEN 2470 

2440 US( I )=UP'T.0/M(1;' 

2450 R< I )=RPC0/MM) 

2452 IF 2A=0 THEN 2460 

2455 PRINT II; I>Mri;KD'' I >>R< I-l ).USa-l > 

2460 GOTO 2610 

2470 IF B=0 THEN LET DS=DT 

2480 IF B=1 THEN LET DS=TS< I >-TS< I-l > 

2490 A=( OM*< P I - AL >*03 )/< 4*< S I NT. < P I - AL >/2 'j.^C-OS< < P I - AL )/2 ) )-< < P I -AL >/2 > ) 
2500 C0=6 . 066* 1 Q'< 22 >*EXP< -8478/KD< I > ) 
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2502 IF 2A=0 THEN 251© 

2505 PRINT 1 1 ; I : MM; KDC I R< I-l >, IJS< I-l ) 

2510 S=EXPC 1500;*:OMA 1 . 38i:*10' < -16 > VR< i-l VKD< I > >-l 
2520 D= . 176*;EXP< -15213 KDC I > > 

2530 LG=LOG< L 'R< I-l >> 

2540 RA=A» C0»S*D/< LGTRC I-l >Z > 

2550 RR=2TRA 

2560 R< I >=R< I-l >+RR:|cR< I-l > 

2570 use I >=US< I-l >+RA4:US<: I-l > 

2575 IF PIJ-0 THEN 2630 
2580 PRINT 

2590 PRINT KD< I ); DS; RA; RR; US< I >; R< I ) 

2600 PRINT A;C0;S;D;LG 
2610 PRINT 

2615 IF PlJ=0 THEN 2630 
2620 PRINT KD<I).L0/I 
2630 NEXT I 

2632 IF II=l»JP THEN 2640 

2635 UP< I H-l , MM >=US< 150 > 

2636 RP< I I+l . MM >=R< 150 ) 

2640 MU<; MM >-< use 150 >-UP< 8 / MM ) >/«JP< 0/ MM > 

2650 MR< MM >=< R< 105 >-RP< 0, MM > >/RP< 0 , MM ) 

2652 SU<MM)=MU<MM> 

2654 SReMM>=MR<MM> 

2655 TU< 1 1 / MM >=MU< MM ) 

2656 TR< I I . MM >=MR< MM > 

2660 NEXT MM 

2665 WAIT 1500 

2700 SU=< SU-e MIN< SU > ) >/< MAX< 3U >-MIN< SU > > 

2710 SR=< SR-< M I N< SR > > >/< MAX< SR >-M IN< SR > ) 

2720 PAGE 

2724 IF II=WP THEN 2730 

2725 IF ZA=0 THEN 2800 

2730 PRINT "MAX D»JTS/UTS= ".MAX<MU) 

2740 PRINT "MIN DUTS/IJTS= 'V.MIH<MU) 


viii-3: 



2750 WINDOW -BH/TH,0. I 
2760 UIEWPORT 250.750.100/600 
2770 SETGR TICS 5.5.GRAT 5.5 
2780 XYPLOT Tr’.SU 
2790 PRINT 
2800 NEXT II 

2809 PAGE 

2810 PRINT "DTS/UTS FOR WP=0.3“ 

2812 PRINT "E''.''D",'‘C'*.'‘B“."A" 

FOR QQ=0 TO WP 

2816 PRINT TlKQQ.0>^TlKQCi!. 1 >.TU<QQ.2).TU<QQ.3>.TU<QQ.4> 
2818 PRINT 
2820 NEXT QQ 

2340 PRINT ”DR/R FOR WP=0.3" 

2842 PRINT "E". "D". "C". "B". "A" 

^46 PRINT*^^< QQ. 0*>. TR< QQ. 1 >. TR< QQ. 2 >. TR< QQ. 3 ). TR< QQ. 4 > 

2843 PRINT 
2350 NEXT QQ 
2855 WAIT 1500 
2860 PRINT 
2865 WAIT 1500 


2900 STOP 

3000 REM COMPUTATION OF HAZ TEMPERATURE HISTORY 

3005 IF 11=0 THEN 3009 

3006 IF MM=LU THEN 3009 

3007 F1=NC 

3008 GOTO 3010 

3009 F1=CC 

3010 P1=P0:»:< 1-F1/100> 

3015 PRINT I I. MM. LU. FI 
3020 P2=P1*< 1-F2/100) 

3030 P1=P1H=>2 

3040 LET N=150 
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3050 DIM XX<N>,TA<N>.TE<N> 

3060 LET TE=LT 

3070 FOR 1=0 TO N 

3080 LET I >=X0+<XB-XA)^NTI 

3090 NEXT I 

3110 FOR 1=0 TO N 

3120 LET X=XX<I> 

3130 GOSUB 3290 
3140 LET TA<I )=ri2 

3145 IF TA< I )>LT THEN LET TA< I )=LT-»-10 
3150 NEXT I 

3160 REM PRINTOUT OF RESULTS 
3170 PAGE 

3175 IF ZA=0 THEN 3280 

3180 PRINT "WELD TQ1PERATURE PROFILE^ Y; •’IJCHES FROM CENTERLINE." 
3190 PRINT "PROFILES TAKEN AT "i2/Wl»100i " *4 PLATE DEPTH." 

3200 FOR 1=0 TO 10 PRINT sNEXT I 
3210 PRINT " TEMP" '.PRINT xPRINT 
3220 PRINT "<DEG F>" 

3230 FOR 1=0 TO 16' PRINT -NEXT I 

3240 PRINT " DISTANCE FROM HEAT SOURCE < INCHES)"; 

3250 WINDOW XA.XB.0,TG 

3260 SETGR WINDOW.TICS 5/ 10/5.5. GRAT 6.6. 3. 3 
3270 XYPLOT XX. TA 
3275 PRINT "-C-W-C-L" 

3280 RETURN 

3290 REM Ktt SUBROUTINE TO COMPUTE TEMPERATURE OF HEAT SOURCE ARRAY 

3300 LET Z2=Z 

3310 LET M2=0 

3320 LET D3=0 

3330 LET D2=2:»cWl 

3340 GOSUB 3650 

3-350 GOSUB 3510 

3360 LET M2=M1+T0+T2 

3370 LET M3=M2 
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3330 LET D3=D3-»'D2 
3390 GOSUB 3430 

3403 IF <M2-M3X.01*M2 THEN GOTO 3420 
3410 GOTO 3370 
3420 RETURN 

3430 REM ttt SUBROUTINE TO ADD NEXT TWO HEAT SOURCES ttt 

3440 LET Z2=D3+2 

3450 GOSUB 3650 

3460 LET M2=M2+Ml 

3470 LET Z2=03-Z 

3480 GOSUB 3650 

3490 LET M2=M2+M1 

3500 RETURN 

3510 REM SUBROUTItC TO COMPUTE LINE HEAT SOURCE TEMPERATURE 

3520 REM DISTRIBUTION 

3530 LET K0=0 

3540 LET T2=0 

3550 LET AR=Ul/2/Al»X 

3560 LET Z3=Ul/2/Al*SG!R<XTX+Y»Y) 

3570 IF AR>88 THEN 3640 

3580 GOSUB 3730 

3590 IF AR<0 THEN 3628 

3600 LET T2=P2/Wl/6.2832/Kl/EXP<AR>*K0 

3610 GOTO 3640 

3620 LET AR=-AR 

3630 LET T2=P2/Wl/6.2332/Kl*EXP<AR>»K0 
3640 RETURN 

3650 REM t*t SUBROUTINE TO COMPUTE TEMPERATURE ttt 
3660 LET Sl=SQR<X4cX+Y*Y+Z2*22) 

3670 IF S1=0 THEN LET M1=UT+1000MF S1=0 THEN 3720 
3680 LET AR=in/2/Al*<Sl+X> 

3690 IF AR>88 THEN LET M1=0 

3700 IF AR>88 THEN GOTO 3720 

3710 LET Ml^Pl/Z/S. 1416/K1/S1/'EXP<AR) 

3720 RETURN 
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3730 REM MODIFIED BESSEL FUNCTION, SECOND KINO, ZEROTH ORDER »»» 

3740 REN m. POLYNOMIAL APPROXIMATION 

3750 IF 23 >2 THEN 3820 

3760 GOSUB 3870 

3770 LET G2=231c23/4 

3780 LET K0=<<7.4E-06:«cG2+l .075E-04)*G2+2.62698E-03>*G2 

3790 LET K0=-;<K0+. 0348859 >*G2+. 2306976 >TG2 

3800 LET K0=<<K0+. 4227842 ):*:G2-. 5772157 >-.5:»:LOG<G2)<cI0 

3810 GOTO 3860 

3320 LET G2=2/Z3 

3330 LET K0=< < 5 . 3208E-04:»cG2~2 . 51 54E-03 >:TG2+5 . 87872E-03 y*GZ 

3840 LET K0=<<K0-. 01062446 >TG2+. 02189568 >T.G2 

3350 LET K0=< <K0-. 07832358 ><:G2+1 253314 >/SQR<Z3VEXP<Z3) 

3860 RETURN 

3870 REM MODIFIED BESSEL RJNCTION, FIRST KIND, ZEROTH ORDER tt* 
3880 REM :♦:*» POLYNOMIAL APPROXIMATION ***■ 

3890 IF 23>3.75 THEN 3950 
3900 LET Gl=23l;23 

3910 LET I0=<<5 923979E-10*G1+6.56017E-08>TGI+6.80123E-06>*G1 
3920 LET I0=<< 10+4. 3394E-04)»G1+. 0156252 >*G1 
3930 LET I0=<I0+.25>*G1+1 
3940 GOTO 3990 

3950 LET I0=< < 153 . 445/23-171 . 822 >/Z3+73 . 2919 VZ3 

3960 LET I0=<a0-15.2595>/23+1.8119S>/Z3 

3970 LET I0=< < 10- . 0830909 >/23+ . 0316855 VZ3 

3980 LET I0=( < 1 0+ . 0498222 >/Z3+ . 3989423 >»EXPC Z3 >/SQR< 23 > 

3990 RETURN 
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